arXiv:astro-ph/0211358vl 15 Nov 2002 


1 


Astrophysical sources of high energy neutrinos 

E. Waxman*^ * 

'^Physics Faculty, Weizmann Institute of Science, Rehovot 76100, Israel 
waxman@wicc.wiezmann.ac.il 

Several high energy, > 100 GeV, neutrino telescopes are currently operating or under construction. Their main 
motivation is the extension of the horizon of neutrino astronomy to cosmological scales. We show that general, 
model independent, arguments imply that ~ 1 Gton detectors are required to detect cosmic high energy neutrino 
sources. Predictions of models of some of the leading candidate sources, gamma-ray bursts and micro-quasars, 
are discussed, and the question of what can be learned from neutrino observations is addressed. 


1. Introduction and summary 


High energy, > 100 GeV, neutrino tele¬ 

scopes are currently operating in deep lake wa¬ 
ter (BAIKAL, http://baikal-neutrino.da.ru/) 
and under Antarctic ice (AMANDA, 
http: / / amanda.berkeley.edu/amanda/). Two 
under-water detectors are currently un¬ 
der construction in the Mediterranean 
(ANTARES, |http://antares.in2p3.fr/] ; NESTOR, 
http:/ /www.cc.uoa.gr^nestor/D , aiming at 
achieving effective volumes ~ O.lkm^, com¬ 
parable to that of AMANDA. Much larger, 
~ Ikm^ telescopes are under construction 
in Antarcti c ice (the IceCube ex tension of 
AMANDA, http://icecube.wisc.edu/ ), and un- 
der development in the Me diterranean (NEMO, 
http: / /nemoweb.Ins.infn.it/ ). For a detailed re¬ 
view of the experiments, and also of their scien¬ 
tific goals, see Q. 

The driving motivation for the construction of 
km-scale neutrino telescopes is the observation of 
cosmic point sources. At present, several neutrino 
telescopes monitor solar MeV neutrinos |Q, and 
may also detect MeV neutrinos from supernova 
explosions, such as supernova 1987A, in our local 
galactic neighborhood. The construction of high- 
energy neutrino telescopes is aimed at extending 
the distances accessible to neutrino astronomy to 
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cosmological scales. This new window onto the 
cosmos will provide a probe of the most powerful 
sources in the universe through observations of 
high-energy neutrinos. 

The existence of extra-Galactic high energy 
neutrino sources is implied by the observations of 
ultra-high energy (UHE), > 10^® eV, cosmic-rays. 
The cosmic-ray spectrum flattens at ^ 10^® eV 
There are indications that the spectral 
change is correlated with a change in composition, 
from heavy to light nuclei . These character¬ 
istics, which are supported by analysis of Fly’s 
Eye, AGASA and HiRes-MIA data, and for which 
some evidence existed in previous experiments [|j , 
suggest that the cosmic ray flux is dominated at 
energies < 10^® eV by a Galactic component of 
heavy nuclei, and at UHE by an extra-Galactic 
source of protons. Also, both the AGASA and 
Fly’s Eye experiments report an enhancement of 
the cosmic-ray flux near the Galactic disk at en¬ 
ergies < 10^®-® eV, but not at higher energies 0. 

Irrespective of the nature of the cosmic-ray 
sources, some fraction of these particles will pro¬ 
duce pions as they escape from the acceleration 
site, either through hadronic collisions with ambi¬ 
ent gas or through interaction with ambient pho¬ 
tons, leading to electron and muon neutrino pro¬ 
duction from the decay of charged pions. 

In § ^ we discuss the upper bound implied 
by UHE cosmic-ray observations on the diffuse 
high energy neutrino flux. In § ^ we show 
that general, model independent arguments im- 
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ply that detectors with masses equal to or larger 
than 1 Gton, equivalent to 1 km^ of water, 
should be constructed in order to detect the ex¬ 
pected neutrino signal in the energy range of 1 
to 10^ TeV. This conclusion applies both for the 
detection of point sources and for the detection 
of diffuse extra-Galactic flux. Even larger detec¬ 
tors may be required at higher energies. The 
construction of detectors with effective volume 
1 km^ at energies » 10^ TeV may require 
the use of techniques different than the under¬ 
water/ice optical Gerenkov technique employed 
by the AMANDA, ANTARES, Baikal, IceCube, 
NEMO, and NESTOR experiments. Such alter¬ 
native technique may be, e.g., the detection of 
coherent radio emission from > 10^® eV neutrino 
induced cascades Q. 

The origin of UHE cosmic-rays is one of the 
most exciting open questions of high energy as¬ 
trophysics The extreme energy of the high¬ 
est energy events poses a challenge to models of 
particle acceleration. Very few known astrophys- 
ical objects have characteristics indicating that 
they may allow acceleration of particles to the 
observed high energies. The detection of high 
energy neutrinos may resolve the puzzle of UHE 
cosmic-ray origin, as the high energy neutrinos, 
unlike the charged cosmic-ray protons which are 
subject to deflection by magnetic fields, will point 
directly to their sources. 

The most powerful known extra-galactic ob¬ 
jects, gamma-ray bursts (GRBs) and active galac¬ 
tic nuclei (AGN), are candidate sources for the 
production of UHE cosmic-rays, and are there¬ 
fore likely sources of neutrinos in the energy range 
of 1 to 10® TeV. GRBs are transient flashes of 
~ 1 MeV gamma-rays lasting typically for 1 to 
100 s, that are observed from sources at cosmolog¬ 
ical distances. The apparent isotropic luminos¬ 
ity of GRBs is ^ 10®^ erg/s. They are believed 
to be powered by the rapid accretion of a frac¬ 
tion of a solar mass of matter onto a newly born 
solar-mass black hole. AGN are persistent sources 
with apparent luminosity reaching ^ 10“^® erg/s. 
They are thought to be powered by mass accre¬ 
tion onto 10®-10® solar-mass black holes, that re¬ 
side at the centers of galaxies. In both GRBs and 
AGN, mass accretion is believed to drive a rela¬ 


tivistic plasma outflow that results in the accel¬ 
eration of high-energy particles, which emit non- 
thermal radiation. A similar process could also 
power Galactic micro-quasars, which may be con¬ 
sidered as a scaled-down version of AGN, powered 
by stellar-mass black holes or neutron stars. 

Despite the overall success of models in ex¬ 
plaining the observed phenomena associated with 
these high energy sources, the models are largely 
phenomenological and our understanding of un¬ 
derlying physical processes is still incomplete. 
In all cases, neutrino observations will provide 
unique information on the physics of the under¬ 
lying engine. 

In § ^ we discuss neutrino emission from GRBs. 
We describe the underlying model and its pre¬ 
dictions for neutrino telescopes, and address the 
question of what can be learned from neutrino 
observations. The discussion illustrates how neu¬ 
trino observations will help resolving open ques¬ 
tions related to the physics that underlies mod¬ 
els of high energy astrophysical sources. In § ^ 
we briefly review neutrino emission from micro¬ 
quasars. 

Finally, it should be emphasized that a detec¬ 
tion of even a handful of neutrino events corre¬ 
lated with GRBs will allow to test for neutrino 
properties, e.g. flavor oscillation and coupling 
to gravity, and to place constraints on deviation 
from Lorentz invariance with accuracy many or¬ 
ders of magnitude better than currently possible. 


We discuss this in some more detail in § 4.6 


2. An upper bound to the diffuse neutrino 
flux 


We first derive in § 2.1 the rate of generation of 
UHE protons implied by cosmic-ray observations. 
We then derive in § 2.2 the upper bound on the 
diffuse neutrino flux. 


2.1. The UHE cosmic-ray generation rate 

Fly’s Eye stereo spectrum is well fitted in the 
energy range ® eV to 10^® ® eV by a sum 
of two power laws: A steeper component, with 
differential number spectrum J oc dom¬ 

inating at lower energy, and a shallower compo¬ 
nent, J oc E“^ ®^, dominating at higher energy. 
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E > 10^® eV. The data are consistent with the 
steeper component being composed of heavy nu¬ 
clei primaries, and the lighter one being composed 
of proton primaries. The observed UHE cosmic- 
ray flux and spectrum may be accounted for by a 
two component, Galactic -I- extra-Galactic model 
|l0[ . For the Galactic component, this model 
adopts the Fly’s Fye fit. 


dn 

dE 


oc F-3-50. 


( 1 ) 


The extra-Galactic proton component is derived 
in this model by assuming that extra-galactic pro¬ 
tons in the energy range of 10^® eV to 10®^ eV are 
produced by cosmologically-distributed sources 
at a rate 


« 3 X lO'^Vg Mpc"^ yr-\ (2) 


with a power law differential energy spectrum 


dup 

dEp 


cx F-" , 


n K, 2. 


( 3 ) 


The corresponding energy per logarithmic decade 
of protons is 


„dhp^ 
P dEp 


0.7 X 10'^'‘erg Mpc"^ yr'^ 


( 4 ) 


The spectral index, n « 2, is that expected for 
acceleration in sub-relativistic collisionless shocks 
in general, and in particular for the GRB model 
discussed in § 4.1. The energy generation rate, 
Eq. ||, is motivated by the GRB energy generation 
rate, compare Eq. ^ with Eq. p^ . 

Figure ^ compares the model prediction with 
the data from the AGASA [Q, Fly’s Eye [||, 
Hires jlTl , and Yakutsk ||l^ cosmic ray experi- 
ments^. The absolute flux measured at 10^® eV 
differs between the various experiments, corre¬ 
sponding to a systematic ~ 11% (~ 19%, ~ 
—7.5%) over-estimate of event energies in the 
AGASA (Yakutsk, HiRes) experiments compared 
to the Fly’s Eye experiment (see also il^,|ir 


■^The Haverah Park data have recently been re-analyzed 
using modern numerical simulations of air-shower develop¬ 
ment [|l3|]. The reanalysis resulted in significant changes of 
inferrea cosmic-ray energies compared to previously pub¬ 
lished results (0 and references quoted therein). This 
improved analysis is available only at energies < 10^® eV. 



Figure 1. Model versus data. The solid curve 
shows the energy spectrum derived from the two- 
component model discussed in Section 2.1. The 
dashed curve shows the extra-Galactic proton 
contribution, Eq. 0 . 


These systematic shifts are well within the sys¬ 
tematic uncertainties quoted by the experiments. 
We have therefore applied the appropriate small 
systematic shifts in absolute energy, to bring the 
various experiments into agreement. The results 
are not sensitive to the choice of absolute energy 
scale [p^ . 

The choice of cosmological model is unimpor¬ 
tant for cosmic ray energies above 10^® eV, which 
is the region of interest here. For definitness, 
a flat universe with = 0.3 and = 0.7, 
and Hubble constant Hq = 65km/sMpc was as¬ 
sumed in Fig. 1^. Also, Eq. ^ represents the local, 
z = 0, energy generation rate. Motivated by the 
evidence for association of GRBs with star for¬ 
mation, it was assumed that the generation rate 
evolves with redshift z like f{z) = (l-bz)^ at low 
redshift, z < 1.9, /(z) = Const, for 1.9 < z < 2.7, 
and an exponential decay at z > 2.7 fl^ . 

The model described above is similar to that 
proposed in |^. The improved constraints on 
UHECR spectrum and flux provided by the re- 
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cent observations of HiRes introduce only a small 
change, compared to |]^, to the inferred energy 
generation rate, given by Eq. ||. 

The suppression of model flux above eV 

is due to energy loss of high energy protons in 
interaction with the microwave background, i.e. 
to the “GZK cutoff” |^. The spectra measured 
by Fly’s Eye, HiRes and Yakutsk are consistent 
with the model, and hence with the existence of 
a GZK cutoff. 

The data from the AG AS A experiment, the ex¬ 
posure of which is ~ 1/3 of the combined expo¬ 
sure of the other experiments, is consistent with 
the model, and with the other experiments, up 
to 10^° eV, but show an excess of events above 
10^° eV. The origin of this discrepancy is not yet 
clear. However, it should be pointed out that 
since the > 10^°eV flux is dominated by sources 
at distances < 50 Mpc, over which the distribu¬ 
tion of known astrophysical systems (e.g. galax¬ 
ies, clusters of galaxies) is inhomogeneous, signifi¬ 
cant deviations from model predictions presented 
in Fig. ^ for a uniform source distribution are ex¬ 
pected at this energy |^. Glustering of cosmic- 
ray sources leads to a standard deviation, cr, in 
the expected number, N, of events above 10^° eV, 
given by a/N = 0.9(do/10Mpc)°'9 |T|, where do 
is the unknown scale length of the source correla¬ 
tion function and do ~ 10 Mpc for field galaxies. 

2.2. The flux bound 

The energy generation rate given by Eq. ^ sets 
an upper limit to the neutrino flux that may be 
produced by sources which, like GRBs and the ob¬ 
served jets of AGN, are optically thin to pion pro¬ 
ducing interactions of protons with source pho¬ 
tons (or ambient nucleons) p^ |. 

If the high-energy protons produced by the 
extra-galactic sources lose a fraction e < 1 of 
their energy through production of pions before 
escaping the source, the resulting present-day en¬ 
ergy density of muon neutrinos is E'^dN^/dEi, « 
0.25etHEpdnp/dEp, where tn « 10^°yr is the 
Hubble time. For energy independent e, the 
neutrino spectrum follows the proton generation 
spectrum, since the fraction of the proton energy 
carried by a neutrino produced through a photo¬ 
meson interaction, E^ 0.05ifp, is independent 


of the proton energy. The 0.25 factor arises be¬ 
cause neutral pions, which do not produce neu¬ 
trinos, are produced with roughly equal probabil¬ 
ity with charged pions, and because in the decay 
of charged pions muon neutrinos carry approxi¬ 
mately half the charged pion energy. Thus, an 
upper limit to the muon neutrino flux [v^ and Dp 
combined) is obtained for e = 1 fl^. 


2^WB 


E < E ^ 

^ — V V 


^■‘^HztH^Ep 


„dh^^ 

Z P 

dEr, 


1.5 X 10"®Cx 


GeV 


cm^ s sr 


( 5 ) 


In the derivation of Eq. || we have neglected the 
redshift energy loss of neutrinos produced at cos¬ 
mic time t < tn, and implicitly assumed that the 
cosmic-ray generation rate per unit (comoving) 
volume is independent of cosmic time. The quan¬ 
tity in Eq. I has been introduce to describe 
corrections due to redshift evolution and energy 
loss. Assuming that the UHE proton, and hence 
neutrino, energy generation rate evolves rapidly 
with redshift, following the evolution of star for¬ 
mation rate [^, see § ^, we find that ~ 3 
(with weak dependence on cosmology). The cor¬ 
rection is small despite the strong evolution with 
redshift, since the universe spends only a small 
fraction of its present age at high z. For no evo¬ 
lution, we have ~ 0.6 < 1 (with weak depen¬ 
dence on cosmology) due to redshift energy loss 
of neutrinos. 

There are two speculative types of sources for 
which the Waxman-Bahcall bound, Eq. H, does 
not apply. These sources, for which we have no 
observational evidence to date, could in principle 
produce a neutrino flux exceeding the Waxman- 
Bahcall limit. The first special type of source 
is one in which neutrinos are produced by pro¬ 
cesses other than photo-meson or proton-nucleon 
interactions. Such sources include, e.g., decays of 
topological defects or of supermassive dark mat¬ 
ter particles (see, e.g., Q). The second type of 
special source is one for which the optical depth 
for photo-production of mesons (or for proton- 
nucleon interaction) is high. Examples of such op¬ 
tically thick scenarios include neutrinos produced 
in the cores of AGNs (rather than in the jets), or 
in a collapsing galactic nucleus [po|. 
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3. Phenomenological considerations: km^- 
scale detectors 


3.1. Detection of point sources 

One may obtain an estimate for the required 
telescope size by considering the minimal flux 
of a source that may be detected by a neutrino 
telescope of effective area A (in the plane per¬ 
pendicular to the source direction) and exposure 
time T. The probability that a muon, produced 
by interaction of a muon neutrino with a nu¬ 
cleon, will cross the detector is given by the ra¬ 
tio of the muon and neutrino mean free paths, 
which (for water, ice) is approximately given by 
= 10-'‘(£;,./100TeV)“, with a=lior < 
100 TeV and a = 0.5 for > 100 TeV. A source 
of energy flux fi, in neutrinos of energy Ei, will 
produce N = {f^/E^)P^^AT events in the detec¬ 
tor. Thus, the flux required for the detection of 
N events is 

u 10-“Af 

X 


E. 


100 TeV 
AT 
km^yr 


erg 


( 6 ) 


For cosmological sources, with characteristic dis¬ 
tance of d ~ c/EIQ ~ 4 Gpc, the minimum lumi¬ 
nosity of a detectable source is therefore 


> 10 ^^ 


E, \ 

100 TeV y 


1 —Q! 


AT 

km^yr J 


erg/s. (7) 


Objects with luminosity > 10"'^^erg/s, more 
than 13 orders of magnitude higher than the Solar 
luminosity, are very rare. The only known class 
of steady sources that produce such high luminos¬ 
ity are AGN. It is therefore clear that km-scale 
neutrino detectors are required for the detection 
of cosmological sources. This argument holds 
also for transient sources: The brightest known 
transients are GRBs, which produce luminosity 
~ 10®^erg/s over ~ 100 s duration. Replacing 
T = 1 yr with T = 100 s in the above equation, 
implies a minimum luminosity Lv ^ 10®^erg/s. 

A lower limit to the source flux is also set by 
the requirement that the signal would exceed the 
background produced by atmospheric neutrinos. 
The atmospheric neutrino flux is approximately 
given by 10-®'5(A,,/500TeV)-^GeV/cm2ssr, 



[GeV] 


Figure 2. The projected sensitivity of km-scale 
detectors, based on Eqs. || and ^ compared with 
the Waxman-Bahcall upper bound on the diffuse 
neutrino intensity, Eq. Also shown are the cur¬ 
rent experimental upper bounds imposed by the 
Baikal and AMANDA experiments and the 
atmospheric neutrino background. 


with f} = 1.7 for E^ < 500 TeV and /3 = 2 
for Ey > 500 TeV. Given that the expected 
angular resolution of the telescopes is approxi¬ 
mately 6 = 1 degree, the source flux for which 
the signal constitutes a 5 (t detection over the 
atmospheric background flux is 


10 “ 


( E, 


- 0.8 


yi0i4.5ev^ 

6 f AT 


Ideg ykm^yr 


- 1/2 


erg 


( 8 ) 


Thus, for km-scale detectors the atmospheric neu¬ 
trino background does not pose more stringent 
constraints on source flux than those imposed by 
the requirement for a detectable signal (except at 
low energies). 


3.2. Detection of diffuse background _ 

Using arguments similar to those used in § 3.1, 
we may obtain the minimum intensity of a diffuse 
neutrino background, that will allow its detection 
by a detector of effective area A observing 27r sr of 
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the sky over a duration T. The intensity required 
for the detection of N events is 




E,, 


100 TeV 
/ AT 
ykm^yr 


l-a 


GeV 
cm^s sr ’ 


(9) 


and the minimum intensity required for a 5cr de¬ 
tection over the atmospheric background is 

El<^, ~ 10 -® 


We 

sitivity of km-scale detectors, based on Eqs. ^ 
and 1^, with the upper bound Eq. || on the dif¬ 
fuse intensity, and with the current experimen¬ 
tal upper bounds imposed by the Baikal and 
AMANDA experiments |2l|. Similar to the case 
of point sources, km-scale detectors are required 
to achieve the sensitivity that may allow detec¬ 
tion of a diffuse background. Note that, as men¬ 
tion in § [^, at energy E^, 10^ TeV the required 

effective volume is ^ Ikm^. 


E,, 


1014.5 ev 
/ AT 


- 0.8 


- 1/2 


GeV 


( 10 ) 


ykm^yry cm^ssr 
re in Figure 0 the projected sen- 


4. GRBs: candidate extra-Galactic 

sources of > 1 TeV neutrinos 

4.1. GRB model and UHE proton produc¬ 
tion 

General phenomenological considerations, 
based on y-ray observations, indicate that, re¬ 
gardless of the nature of the underlying sources, 
GRBs are produced by the dissipation of the 
kinetic energy of a relativistic expanding fireball 
A compact source, tq ~ 10^ cm, produces 
a wind, characterized by an average luminosity 
L ~ lO^^erg s“^ and mass loss rate M. At small 
radius, the wind bulk Lorentz factor, T, grows lin¬ 
early with radius, until most of the wind energy 
is converted to kinetic energy and T saturates at 
r ~ LjMc? ~ 300. Variability of the source on a 
time scale At ~ 10 ms, resulting in fluctuations 
in the wind bulk Lorentz factor T on a similar 
time scale, results in internal shocks in the ejecta 
at a radius r ~ « E^cAt ^ tq. It is assumed 


that internal shocks reconvert a substantial part 
of the kinetic energy to internal energy, which 
is then radiated as y-rays by synchrotron and 
inverse-Compton radiation of shock-accelerated 
electrons. At a later stage, the shock wave driven 
into the surrounding medium by the expanding 
fireball ejecta leads to the emission of lower- 
energy ’’afterglow.” 

Since the observed radiation is produced at a 
large distance, 3> tq, from the underlying source, 
the model is largely independent of the nature 
of the underlying compact object, as long as this 
object is capable of producing a wind with the 
properties implied by observations. The underly¬ 
ing GRB progenitors are yet unknown. Collapse 
of massive stars and mergers of compact objects 
(e.g. of binary neutron stars) are the most widely 
discussed candidates 

The observed radiation is produced, both dur¬ 
ing the GRB and the afterglow, by synchrotron 
emission of shock accelerated electrons. In the re¬ 
gion where electrons are accelerated, protons are 
also expected to be shock accelerated. The inter¬ 
nal shocks within the expanding wind are mildly 
relativistic (in the wind rest frame), and hence 
we expect our understanding of non-relativistic 
shock acceleration to apply to the acceleration of 
protons in these shocks. In particular, the pre¬ 
dicted energy distribution of accelerated protons 
is expected to be dup/dEp oc similar to 

the predicted electron energy spectrum, which is 
consistent with the observed photon spectrum. 

Several constraints must be satisfied by wind 
parameters in order to allow proton acceleration 
to high energy Ep. We summarize below these 
constraints. The reader is referred to [E4 25 


for a more detailed discussion. The requirement 
that the acceleration time be smaller than the 
wind expansion time sets a lower limit to the 
strength of the wind magnetic field, which may 
be expressed as a lower limit to the ratio of mag¬ 
netic field to electron energy density H, 


ub/us > O.O2T2 ^Ep 2qL^^^2j 


( 11 ) 


where Ep = 10^^Ep^20 eV, T = 10^-®r2.5 and 
Lj = 10 ®^L.y_ 52 erg/s is the wind y-ray luminos¬ 
ity. A second constraint is imposed by the re¬ 
quirement that the proton acceleration time be 
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smaller than the proton energy loss time, which is 
dominated by synchrotron emission. This sets an 
upper limit to the magnetic field strength, which 
in turn sets a lower limit to T 




( 12 ) 


Here, At = 10“^At_2 s. As explained in p4[ , 
the constraints Eq. and Eq. ^ hold regardless 
of whether the fireball is a sphere or a narrow 
jet (as long as the jet opening angle is > l/E). 
The luminosity in Eq. is the ’’isotropic equiv¬ 
alent luminosity”, i.e. the luminosity under the 
assumption of isotropic emission. 

Internal shocks within the wind occur over a 
wide range of radii, corresponding to a wide range 
of variability time scales: At ~ 1 ms, the dynam¬ 
ical time of the source, to At ~ 1 s, the wind 
duration. Protons are therefore accelerated to 
high energy over a wide range of radii. In par¬ 
ticular, at large radii the external medium affects 
fireball evolution, and a ’’reverse shock” is driven 
backward into the fireball ejecta and decelerates 
it, due to the interaction with the surrounding 
medium. This shock is also mildly relativistic, 
and its parameters are similar to those of an inter¬ 
nal shock with At ~ 10 s. Protons may therefore 
be accelerated to ultra-high energy in this shock 
as well 

and Eq. 


The constraints Eq. 
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ably similar to those inferred from y-ray obser¬ 
vations: r > 300 is implied by the y-ray spec¬ 
trum in order to avoid high pair-production op¬ 
tical depth, magnetic field close to equipartition, 
us/ue ~ 0.1, is required in order to account for 
both y-ray emission and afterglow observations 
1^ . This suggests that GRBs and UHE s may 
originate from common sources. 

The suggested association between GRBs and 
the sources of UHE protons is further strength¬ 
ened by comparing the proton energy generation 
rate, Eq. with the GRB energy generation rate. 
The evidence for association of GRB sources with 
star-formation |Q, suggests that the GRB rate 
evolves with redshift following the star-formation 
rate. Under this assumption, the local (z = 0) 
GRB rate is w 5 x 10“^°Mpc“^ yr“^ and 
their average 0.1 MeV to 2 MeV y-ray energy re¬ 
lease is 2.5 X 10®^ erg [M, corresponding to an 


energy generation rate ol|| 

« 1.3 X lO^^erg Mpc“3 yr-b (13) 

This rate is remarkably similar to the energy gen¬ 
eration rate of UHE protons inferred from cosmic- 
ray observations, Eq. ||. Note, that although the 
proton generation rate, Eq. ||, is approximately 
twice the GRB y-ray generation rate, Eq. |3[ the 
corresponding energy per logarithmic decade of 
electrons implied by Eq. 

jaGRB 

~ lO'^^erg Mpc-3 yr-b (14) 

ClSq 

is similar to Eq. Thus, GRBs may be the 
sources of observed UHE cosmic-rays, provided 
they produce similar energy in MeV y-rays, or, 
equivalently in high energy electrons, and in high 
energy protons. 

4.2. ~ 100 TeV neutrinos 

Protons accelerated in the fireball to high en¬ 
ergy lose energy through photo-meson interaction 
with fireball photons. The decay of charged pions 
produced in this interaction results in the produc¬ 
tion of high energy neutrinos. The key relation is 
between the observed photon energy, , and the 
accelerated proton’s energy, Up, at the threshold 
of the A-resonance. In the observer frame, 

E^ Ep = 0.2 GeV^ T^ . (15) 

For r Ri 300 and E^ = 1 MeV, we see that charac¬ 
teristic proton energies ^ 10^® eV are required to 
produce pions. Since neutrinos produced by pion 
decay typically carry 5% of the proton energy, 
production of ~ 10^^ eV neutrinos is expected 



The fraction of energy lost by protons to pions, 

u is ©HI] 

Mep) ^ 0.2L\%At Z2^^. (16) 

^The quoted GRB rate density and 7 -ray energy assume 
that the GRB emission is isotropic. If emission is confined 
to a solid angle AQ < 47 r, then the GRB rate is increased 
by a factor {Ar 2 / 47 r)~^ and the GRB energy is decreased 
by the same factor. However, their product, the energy 
generation rate, is independent of the solid angle of emis¬ 
sion. It is currently believed that GRH’s are beamed on 
average into a solid angle of 47r/500 [M, which implies a 
rate higher (total energy lower) by a factor of 500 com¬ 
pared to values inferred assuming isotropic emission. 





Assuming that GRBs generate similar energy 
in high-energy protons and electrons, account¬ 
ing therefore for the observed UHE cosmic-ray 
flux, then using Eq. ||, the expected GRB muon- 
neutrino flux is 

« 0.9 X 10-®^GeVcm-2s-isr-b(17) 

This neutrino spectrum extends to ~ 10^® eV, 
and suppressed at higher energy due to energy 
loss of pions and muons. Comparing Eq. with 
Eq. 1^, we find that ~ 20 neutrino-induced muon 
events per year are expected (over 47r sr) in a 
cubic-km detector. Note, that GRB neutrino 
events are correlated both in time and in direc¬ 
tion with gamma-rays, and hence there detection 
is practically background free . 

4.3. ~ 10^^ eV ’’Afterglow” neutrinos 

High energy neutrino emission may also result 
from photo-meson interactions of protons acceler¬ 
ated to high energies in the reverse shocks driven 
into the fireball ejecta at the initial stage of in¬ 
teraction of the fireball with its surrounding gas, 
which occurs on time scale T ~ 10 s (see § E3)- 
Optical-UV photons are radiated by electrons ac¬ 
celerated in shocks propagating backward into 
the ejecta. The interaction of these low energy, 
10 eV-1 keV, photons and high energy protons 
produces a burst of duration T of ultra-high 
energy, lO^'^-lO^® eV, neutrinos, as indicated by 

Eq.im. 

The expected muon neutrino intensity depends 
on the density of the surrounding medium [^7U32| , 
which may differ by orders of magnitude between 
models assuming different GRB progenitors. If 
GRB fireballs expand into typical density inter¬ 
stellar medium, n ^ lcm“^, as expected if the 
GRBs are the result of mergers of compact ob¬ 
jects, the expected intensity is 

« 10-1° (^^^^^^GeVcm-^s-isr-i, (18) 

with 13 = 1/2 for > lOU’eV and (3 = 1 other¬ 
wise. If GRB fireballs expand into massive stellar 
winds (n ^ 10^cm-°), as expected if GRBs result 


from the collapse of massive stars, then 

El<l>, « 10-Smin{l,£;,/I0i^eV}-^^. (19) 

cm“^ssr 

The neutrino flux is expected to be strongly sup¬ 
pressed at energy > IQi® eV, since protons 
are not expected to be accelerated to energy 
Ep > 10^0 eV. 

4.4. ^ 5 TeV ’’Collapsar” neutrinos 

Lower energy, ~ 5 TeV, neutrino emission may 
be expected in the case where GRBs originate 
from the collapse of massive stars [Q. In this 
scenario, accretion onto a black hole, produced 
by the collapse of the massive star (the ’’collap¬ 
sar”) core, drives a relativistic jet that propagates 
through the stellar envelope along the collapsar 
rotation axis. The shocks producing the q-rays 
must occur after the fireball has emerged from 
the stellar envelope. While the jet is making its 
way out of the star, its rate of advance is slowed 
down in a termination shock that heats the stellar 
plasma to keV temperatures, and additional in¬ 
ternal shocks are expected in the pre-deceleration 
interior jet. The latter can accelerate protons to 
> 10° GeV, which interact with the X-ray pho¬ 
tons in the stellar jet cavity leading to pion pro¬ 
duction and hence electron and muon neutrinos 
(and anti neutrinos) with energies E^, > 5 TeV. 
These neutrinos appear as a precursor signal, last¬ 
ing for time scales of tens of seconds prior to the 
observation of any q-rays produced outside the 
star by a collapsar-induced GRB. The TeV neu¬ 
trino fluence from an individual collapse at cos¬ 
mological distance z ~ 1 implies ~ 0.1 upward 
moving muons per collapse in a 1 km° detector. 


4.5. Implications: GRB models 

The emis sion of ^ 100 TeV neutrinos dis¬ 
cussed in § |4.2| is independent of the underly¬ 
ing GRB progenitor. It is a natural outcom e of 
the ’’generic” fireball model described in § 4.1, 
where the observed q-rays are produced by inter¬ 
nal shocks within a dissipative relativistic wind. 
The major underlying assumption, upon which 
the predictions depend, is that the fireball mo¬ 
mentum is carried (after the initial stage of ac¬ 
celeration) by relativistic protons. The predicted 
flux, Eq. 1^, implies ^ 20 neutrino-induced muon 
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events per year in a cubic-km detector. Neutrino 
telescopes may therefore allow to test underlying 
assumptions of the fireball model. Detection of 
the predicted signal will also provide strong sup¬ 
port for the model of UHE cosmic-ray production 
in GRBs. 

The emission of ~ 10^^ eV neutrinos, § 4.3 


is also a natural prediction of the fireball 
model. The expected neutrino intensity depends 
strongly, however, on the GRB progenitor, com¬ 
pare Eqs. 1^ and 19. Detection of such UHE neu¬ 
trinos may therefore provide constraints on the 
progenitor type. 

Finally, emission of ~ 5 TeV ’’collapsar” neu¬ 
trinos is expected only for collapsar GRB pro¬ 
genitors. The detection of a ~ 5 TeV neutrino 
precursor will therefore provide a clear signature 
of the collapsar model. 


4.6. Implications: Lorentz invariance, The 
-weak equivalence principle, and neu¬ 
trino oscillations 

Detection of neutrinos from GRBs could be 
used to test the simultaneity of neutrino and pho¬ 
ton arrival to an accuracy of ~ 1 s (~ 1 ms for 
short bursts), checking the assumption of spe¬ 
cial relativity that photons and neutrinos have 
the same limiting speed [^. These observations 
would also test the weak equivalence principle, 
according to which photons and neutrinos should 
suffer the same time delay as they pass through 
a gravitational potential. With 1 s accuracy, a 
burst at a distance of 1 Gpc would reveal a frac¬ 
tional difference in limiting speed ~ 10“^^, and 
a fractional difference in gravitational time de¬ 
lay of order 10“® (considering the Galactic poten¬ 
tial alone). Previous applications of these ideas 
to supernova 1987A, where simultaneity could 
be checked only to an accuracy of order several 
hours, yielded much weaker upper limits: of or¬ 
der 10“® and 10“^ for fractional differences in the 
limiting speed and time delay, respectively. 

The model discussed above predicts the pro¬ 
duction of high energy muon and electron neu¬ 
trinos with a 2:1 ratio (tau neutrinos may 
be produced by photo-production of charmed 
mesons; however, the higher energy threshold 
and lower cross-section for charmed-meson pro¬ 


duction, compared to pion production, typically 
imply that the ratio of charmed-meson to pion 
production is 10“'* |Q). Because of neutrino os¬ 
cillations, neutrinos that get here are expected to 
be almost equally distributed between flavors for 
which the mixing is strong. Upgoing taus, rather 
than muons, would be a distinctive signature of 
such oscillations. It may be possible to distin¬ 
guish between taus and muons in a km^-scale de¬ 
tector, since at 10^ TeV the tau decay length is 

1 km. This will allow a ”tau appearance ex¬ 
periment” . 

5. Micro-quasars: Galactic candidate 

sources 

The jets associated with Galactic micro¬ 
quasars Q are believed to be ejected by ac¬ 
creting stellar-mass black holes or neutron stars. 
Much like for AGN, the content of the jets is an 
open issue. The dominant energy carrier in the 
jet is at present unknown (with the exception of 
the jet in SS433). Scenarios whereby energy ex¬ 
traction is associated with spin-down of a Kerr 
(rotating) black hole favor electron-positron com¬ 
position (although baryon admixture is an issue), 
while scenarios in which an initial rise of the X-ray 
flux leads to ejection of the inner part of the accre¬ 
tion disk imply electron-positron jets, as widely 
claimed to be suggested by the anti-correlation 
between the X-ray and radio flares seen during 
major ejection events. A possible diagnostic of 
electron-positron jets is the presence of Doppler- 
shifted spectral lines, such as the Ho lines as seen 
in SS433. The detection of such lines from jets 
having a Lorentz-T factor well in excess of unity 
(as is the case in the super-luminal micro-quasars) 
may, however, be far more difficult than in SS433, 
as the lines are anticipated to be very broad. 

Neutrino telescopes may prove to be a much 
sharper diagnostic tool. If the energy content of 
the jets in the transient sources is dominated by 
electron-proton plasma, then a several hour out¬ 
burst of 1 to 100 TeV neutrinos (and high-energy 
photons) produced by photo-production of pions 
should precede the radio flares associated with 
major ejection events Several neutrinos may 
be detected during a single outburst by a km^- 
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scale detector , thereby providing a power¬ 

ful probe of micro-quasar jet physics and of their 
innermost structure. 
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